Secondary Coulomb Blockade Gap in a Four-Island Tunnel-Junction Array 
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In the ring-shaped tunnel-junction array with four islands, the secondary Coulomb blockade gap 
in a low bias-voltage range is observed in the I — V characteristics. We attribute its appearance to 
the unique topology of the array which induces up to two electrons to get trapped inside. We have 
analyzed the formation and destruction of the gap in terms of detailed single-electron tunneling 
processes. The negative differential resistance behavior when the thermal and quantum fluctuations 
are present is also studied. 
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I. INTRODUCTION 

In a linear one-dimensionaL CAD) tunnel- j unction ar- 
rays with identical junctions,uT3 electrons have no dif- 
ficulty in traveling through the array once they are in- 
jected into the arrayB In other words, electrons do not 
get trapped inside the array, unless the characteristics 
of some junctions of the array (e. g. junction resistance 
or capacitance) are artificially altered so as to trap the 
electrons. On the other hand, the ring-shaped array that 
is composed of two ID arrays between the source and 
the drain electrodes provides a natural way to trap the 
electrons inside the array. In the four-island array of 
Fig. |], which has the smallest possible size for the ring- 
shaped array, one electron may get trapped on island 1 
(the island that is closest to the source electrode) or two 
electrons, on islands 2 and 4 (the top and the bottom 
islands) .0 We have recently reported that up to six elec- 
trons .can be trapped inside the sixty-island ring-shaped 
array.ta 

An important consequence of electrons' getting 
trapped inside the ring-shaped array is appearance of 
the multiple Coulomb blockade (MCB) gaps in the 
zero-temperature I — V characteristics.! While electrons 
are trapped inside, they block the electrical conduction 
through the array so no current flows. Let us denote such 
zero-current region with n t trapped electrons by 72 (n t ), 
and let us suppose that up to Nt electrons can be trapped 
inside. That multiple electrons may get trapped implies 
that there exist multiple regions 72.(1), 72.(2), • • •, !Z(N t ) 
in the I —V curve where no current flows. If all of the in- 
sulating regions are connected adjoiningly, a single gap in 
the form of an extended Coulomb blockade (CB) gap will 
be seen in the I—V characteristics. (Recall that a normal 
CB gap appears in the region where there is no electron 
inside the array: i. e. 72(0) by our notation.) On the 
other hand, if neighboring regions !Z(n t ) and 72(rt t + 1) 
are separated by a conducting region, the MCB gaps will 
be seen in the I— V characteristics. We have shown that, 
in the ring-shaped array, during the transition from odd 
to even n t there may exist such a conducting region, but 
during the transition from even to odd n t no such con- 



ducting region exists E We qualitatively explained the be- 
havior by noting that the topology of the array demands 
that dynamically unstable charge configurations must be 
passed during the odd-to-even transition, but not neces- 
sarily during the even-to-odd transition. Our argument 
leads us to state that the MCB gaps are a unique trans- 
port property of the ring-shaped array, which cannot be 
seen in linear ID arrays, even with artificially trapped 
electrons. 

In this paper, we will quantitatively describe the 
physics of the single-electron tunneling in the ring-shaped 
array, by taking the four-island array as the prototype 
system. Due to its simplicity, we will be able to identify 
the detailed tunneling processes which leads to the for- 
mation and destruction of the secondary Coulomb block- 
ade (SCB) gap (note that for the four-island array, only 
one extra CB gap is possible, which we denote by the 
SCB gap.) Then we will be able to discuss under which 
conditions the SCB gap appears, determine the precise 
location of the gap, and predict the evolution of the gap 
as the thermal and quantum fluctuations are introduced. 
Based on the discussions, the physical properties of the 
MCB gaps of larger size array may be similarly under- 
stood. We will also investigate the effect of the uniform 
background charges on each island on the SCB gap. 



II. THEORY 

Our array system consists of four metallic or semicon- 
ducting islands, positioned as depicted in Fig. [I], between 
the source and the drain electrodes. Nearest-neighboring 
islands are connected by tunnel junctions so that there 
are two conducting paths between the source and the 
drain: upper one via the top island (island 2) and lower 
one via the bottom island (island 4) . Let us assign C and 
R for the tunnel-junction capacitance and resistance, re- 
spectively, and C x for the cross-island capacitance (be- 
tween islands 1 and 3 and between islands 2 and 4). 
Furthermore, each island is capacitatively coupled to the 
ground with same capacitance Co- The capacitance ma- 
trix A is then written as 
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r± = (17 ± Vl7)/2. 



(7) 




FIG. 1. The four-island array. The tunnel-junctions are 
indicated by dotted lines. 



The current / is calculated for the constant applied 
voltage V by using the standard Monte Carlo (MC) 
method witk n the transition rates determined by the 
Golden-rule;aQ the transition rate T for a single-electron 
tunneling over a junction is given by 
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-AF 



Re 2 l-exp(AF/fc B T) : 



(8) 



where AF is the free energy change by the tunneling 
event. The free energy consists of the electrostatic energy 
contribution and the work done by the external voltage 
source as usual. One may conveniently obtain the free en- 
ergy change AF by a tunneling process involving islands 
i and j through the relationshipo 



AF 



(9) 



A = 
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(1) 



where D 1 = -3C - C x - C and D = -2C - C x - C . 

When a constant voltage V is applied between the 
source and the drain and n^e charges are there on is- 
land i (rii is an integer), the potential 4>i of island i is 
given by, in the matrix form, 

( /» = A- 1 Q, (2) 

where the column vector Q is defined by 

Qi = qi - CVS itl , (i=l,..,4) (3) 

where = rc^e. By inverting the capacitance matrix 
given by Eq. (pi), we have 




(4) 



where 



where 



a± = A+/t_ + A_/t+ ± Ai/2, 
6± = A+/t + + A_/t_±A 2 /2, 
c= (A_ -A+)/2Vl7, 



Ar 1 =D'- C x , 
A 2 1 = D — C x , 

D + C x -(l± VT7)/2, 



(5) 



(6) 



A ± X 



and 



where (pi and <p\ are potentials of «th island before and 
after the tunneling, respectively. Since only one electron 
is involved in a single tunneling process, the charges on 
each island are changed by ±e at most, so we can easily 
calculate AF via Eqs. @)-@) and Eq. (g). 

We also used the master equation approach assisted by 
the MC method as follows, especially when tace events 
should be taken care of. The master equationolil] is writ- 
ten as 



rlP N * N 



j i 1 



(10) 



where Tjj is the transition rate from state j to state i 
and Pi is the probability of finding state i among total 
N* states considered in the calculation. The states con- 
sidered in Eq. ( |To| ) are composed of frequent states that 
are readily obtained by a MC run and of rare states, if 
necessary, which may have to be inserted manually with 
some physical reasoning if they are hardly visited by the 
MC runs. Eq. (|l0|) is solved either iteratively or by di- 
rectly solving the equilibrium-state equation: 



GP = 



(11) 



matrix where Gij = Ty 



where G is a N* by N* 

Sij J2i Tii an d P is the vector of the state probabilities. 
Since the matrix G is singular, we take a sub-matrix H 
of size N* — 1 by removing row k and column k from the 
full matrix G. Then 



/ Pi \ 

Pk-i 
Pk+i 

V jv J 



F fc H- 



Gk-l,k 



(12) 
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where Pk is obtained by 



N*-l 

V = l- E H ij lG 3+k,k- (13) 

In the master equation approach, a priori knowledge 
about the relevant states is the key to the successful cal- 
culation. ■— I 

We have also considered cotunnelingjlil up to the sec- 
ond order. The reason that only up to the second-order 
is considered here is that we will be interested in cotun- 
neling out of the stationary charge configuration, as will 
be discussed later, where the charges are trapped on the 
top (island 2) and bottom (island 4) islands of the array. 
The cotunneling rate is then given by the exact double- 
junction formulaO 

= Rk l2-(l+ 2 AFjAF 2 \ 
7 4n 2 e 2 R 1 R 2 \ \ AF AF — AFi — AF 2 ) 

x f E ln ( : ~ Ai 7 AF j | A ^ ( 14 ) 

where AFi and AF 2 are the free energy changes for the 
intermediate tunneling processes which are involved in 
the cotunneling process, and AF is that for the cotun- 
neling. In the above equation, Rk — h/e 2 ~ 25.8KO, 
and Ri and R 2 are the resistances of the junctions in- 
volved in the cotunneling process. 

In this paper, the units of the current, voltage, and 
temperature are I = e/RC, V — e/C, and T — e 2 /kBC, 
respectively. 

III. RESULTS AND DISCUSSIONS 

A unique feature in the I—V characteristics of the four- 
island tunnel-j unction array is the existence of another 
CB gap in a low voltage region at low temperatures,Q 
as shown in Fig. |[. In the figure, where Cq/C = 3, 
C x /C = 0, and T/T = 0, the SCB gap in the interval of 
[0.16,0.20] is seen in addition to the primary CB gap in 
the interval of [0, 0.13]. The origin of the additional gap is 
attributed to the formation of the trapped-electron con- 
figuration or the stationary charge configuration (SCC) 
where two electrons trapped on the topjj-nd the bottom 
islands block the electrical conduction,!™ as schemati- 
cally shown in inset of Fig. ||. 

A SCC is a local minimum of the free energy in the 
configuration space. That is, for each transition from 
the SCC to its adjacent configurations (the configurations 
which can be reached by a single-electron tunneling from 
the SCC) the free energy change AF > 0; otherwise, the 
electrons will not be trapped. In the following we will 
describe the properties of the SCCs in the four-island 
array. 
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FIG. 2. The I — V characteristics of the four-island array 
for Co/C = 3 and C x /C — at zero temperature (thick solid 
line). In the inset is shown the stationary charge configuration 
responsible for the secondary Coulomb blockade gap in the 
interval [0.16,0.20]. The units of the current and the voltage 
are 7/7 and V/V, respectively. 



A. Stationary Charge Configurations 

There are two SCCs that are noteworthy in the four- 
island array. The first one is the SCC with one trapped 
electron on island 1 (SCC-1). If we denote a charge 
configuration by { — ti^} where m is the number of elec- 
trons on island i, SCC-1 corresponds to {—1, 0, 0, 0}. One 
of distinct features of the ring-shaped four-island array 
in contrast to the linear ID arrays is the existence of 
the SCC-1. In linear ID arrays of four islands with 
identical junctions, the tunneling process which triggers 
the conduction is {0,0,0,0} — * { — 1,0,0,0}; once an 
electron tunnels onto the first island connected to the 
source, it has no problem in tunneling through the ar- 
ray. More specifically, let us denote V\ as the thresh- 
old voltage for entrance of one electron onto the array 
({0, 0, 0, 0} {-1, 0, 0, 0}; let us denote the free energy 
change of the process by AFi) and Vj as the threshold 
voltage for the electron to move forward ({—1, 0, 0, 0} — > 
{0,-1,0,0,}). In linear ID arrays, Vt < Vi, which im- 
plies that at the threshold voltage Vi, the free energy 
change AF t = F{0, -1,0,0} - F{-1, 0,0,0} is already 
less than zero, so the threshold voltage for electrical con- 
duction is V\ in that case. In the ring-shaped four-island 
array, however, V% < Vt, as will be seen shortly, which 
implies |that there is a range [Vi,Vt] where the SCC-1 
persistsJla 

AFi and AF t of our four-island array are given by, 
using Eqs. (§)-(§) and Eq. (g), 

AFi = F{-1, 0, 0, 0} - F{0, 0, 0, 0} 

= -(l + o + )V-o+/2 (15) 

and 
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AF t = F{0, 



-1, 0,0} -F{-1,0, 0,0} 
a + )V - (b + - o+)/2. 



(16) 



The threshold voltages Vi and Vj are zeros of the above 
free energy differences. For C x /C = (the strong screen- 
ing case), for simplicity, we get from Eqs. (pjl) and (fill) 



and 



V 



(C + l)(C +4) 
2(C (C + 3)(C + 4) + 2) 



(Cq + 1)(Cq+4) 
2(Co + 2)(C 2 +4Co + ir 



(17) 



(18) 



From Eqs. ( |l7j ) and (|l8|), one can easily see that V t > V\ 
for all Co > 0. The same is true for arbitrary C x /Cq. We 
thereby show that SCC-1 exists in the bias- voltage range 
of [V\, Vt] in our four- island array. 

The second SCC, which is at the heart of our discus- 
sion here, is the one with two trapped electrons on is- 
lands 2 and 4 (SCC-2), whose charge configuration is 
{0, —1, 0, —1}. To investigate the formation and destruc- 
tion of the SCC-2, we need to consider following two tun- 
neling processes: 1) the tunneling process {0, —1, 0, 0} — > 
{ — 1,-1,0,0} which leads to onset of the SCC-2 (once 
the configuration { — 1,-1,0,0} is reached, the desired 
stationary configuration will be eventually reached from 
there because the tunneling process { — 1,-1,0,0} — > 
{0, —1, 0,-1} is a frequent process) and 2) the tunneling 
process {0,-1,0,-1} — * {0,-1,-1,0} (or, equivalently, 
{0, -1, 0,-1} -> {0, 0,-1, -1}) which leads to break-up 
of the SCC-2. Let us denote the free energy changes of 
the former and the latter tunneling processes as AF 2 and 
AFb, respectively. Then, 



AF 2 = F{— 1, — 1, 0, 0} — F{0, —1, 0, 0} 
= -(l + o+)V- (a+ + 2c)/2 



(19) 



and 



AF b = F{0, —1,-1,0} — F{0, — 1, 0, —1} 

= (c - a_)V + ((o+ - b+) + 2(c - b-))/2. (20) 

At the threshold voltage V 2 where AF 2 of Eq. (fij) 
becomes zero, a second electron can now enter the ar- 
ray, and after subsequent tunneling events, the charge 
configuration {0,-1,0,-1} will be eventually reached. 
From the configuration, the only possible tunneling pro- 
cess is {0,-1,0,-1} — > {0,-1,-1,0} (or equivalently, 
{0,-1,0,-1} -> {0,0,-1,-1}). If the free energy 
change for the process AFb is greater than zero, the 
tunneling process is energetically unfavorable so that the 
electrons get trapped and consequently the SCC-2 is es- 
tablished. See Fig. ||. Thus the necessary condition for 
the formation of the SCC-2 is V% < Vb, where Vb is the 
zero of AFb ■ 



Voltage 




FIG. 3. The free energy differences AF t , AF 2 , and AFb 
(thin solid lines) are shown together with the I — V curve of 
Fig. ^| (thick solid line, left axis). Also shown is the average 
number of electrons (n) (dotted line, right axis). The thresh- 
old voltages Vi, Vt, V2, and Vb are marked. The unit of the 
free energy differences is e 2 /C. 



There are two tunneling processes which lead to de- 
struction of the SCC-2. The first one is the tunneling pro- 
cess where the trapped electrons are forced to move for- 
ward: {0, —1,0, —1} — * {0, —1, —1, 0} which is energeti- 
cally favorable if AF b < 0, that is, if V > V b . The other 
one is the tunneling process which introduces a third elec- 
tron into the array: {0,-1,0,-1} — > { — 1,-1,0,-1}, 
which is energetically possible only when V > V3 , where 
V3 is the threshold voltage for the process. Therefore, 
the SCC-2 persists in the voltage range [V2, min(Vh, V3)]. 

Therefore, the SCC-2 is established if the following 
condition is satisfied: 



Vt<V 2 <mm(V b ,V 3 ). 



(21) 



Then a current peak in the voltage interval of [Vt,I^] 
and the SCB gap in the interval of [V2, min(Vb, V3)] will 
be seen in the I — V at zero temperature. 

Let us now discuss when the condition of Eq. j2l| ) 
holds. For C x /C = 0, for simplicity, threshold voltages 
V2, Vb and V3 are given by 



V 2 = 



(Co + 3)(C + 4)-2 

2(C (O)+3)(Cu + 4) + 2) : 
C% + Co - 4 



and 



2(C + 2)(C + 1)' 



V 3 « V 2 + (C/C eff ) 2 e/C 



(22) 
(23) 

(24) 



where C e// = ^Cl + AC C. From Eqs. (fij) and (||)- 
(0) , we obtain that the equality of Eq. (|2l|) is satisfied 
if Cq > Co ~ 2.7. In general, when the cross capacitance 
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C x is nonzero, the expressions for the threshold voltages 
are much more complex than in the strong screening case, 
and one has to resort to numerical calculation for the 
evaluation of the equality of Eq. (^l]). In Fig. |[ we show 
Cq, the critical value of Cq for the observation of the 
secondary Coulomb blockade gap, with respect to the 
screening factor C x . The figure shows that as the cross 
capacitance C x increases Co also increases. 

In summary, we have shown that the SCC-1 exists 
in the range 72.(1) = [Vi,V t ] and the SCC-2 may exist 
in the range 72.(2) = [V 2 , min(Vf,, V3)], if the condition 
Vt < Vi < mixx(Vb, V3) is satisfied. The regions 72.(1) and 
72.(2) are zero-current regions, but in between the two re- 
gions, the current flows such that it appears in the I — V 
characteristics that the SCB gap shows up in the range 
72.(2). In terms of number of trapped electrons n t , we 
observe that the system becomes conducting during the 
transition from nt = 1 to 2 (see Fig. ^ where the average 
number of electrons inside the array (n) is also shown). 
Noting that both the SCC-1 and the SCC-2 are stable 
charge configurations in time, we remark that the reason 
that the system becomes conducting during the transi- 
tion is that dynamically unstable charge configurations 
in time must be passed en route. If we start with the 
SCC-1 and try to reach the SCC-2, the charges should 
physically move in the meantime and they keep tunnel- 
ing through the array until the charges are balanced as 
required by the geometry of the array (or electrostatic 
force). In short, the SCB gap results because the topol- 
ogy of our ring-shaped array makes it possible that unsta- 
ble conducting state exists between two stable insulating 
states. For the array of larger size having the same topol- 
ogy (i.e. the array with two branches between the source 
and the drain), this behavior is further generalized: that 
is, the current peaks can appear during each transition 
from odd (n) to even (n)B We have reported in Ref. || 
that three peaks can be observed for the array of 60 is- 
lands, during each transition of (n) = 1^2,3^4, and 
5^6. The reason that the current peaks arise only when 
nt changes from odd to even is that the transition from 
even to odd nt does not involve unstable charge states 
during the transition. 



B. Effect of Uniform Background Charges 

We have so far dealt with the system with neutral back- 
ground charge on each island: i. e. qi — n^e in Eq. (|^). 
A uniform background charge qo on each island may be 
induced by attaching a (metallic) plate close to the array 
and applying a voltage V g = qo/Co between the plate 
and the ground. Then qi = r^e + qo in Eq. (||), 

The effect of the uniform background charge qo on the 
threshold voltages introduced above is to shift them by a 
certain amount linearly proportional to qo but at different 
rates for different threshold voltages. For example, the 
threshold voltages V 2 and VJ, become: 




c/c 



FIG. 4. The critical value Cq/C for the observation of the 
SCB gap at zero temperature with respect to the cross capac- 
itance C x /C for qo/e = 0.1 (crosses), (solid line), -0.1 (solid 
circles), and -0.2 (open circles), qo is the uni form background 
charges induced on each island (see Section [II Bj ). 



V 2 {q ) =V 2 (q =0) + q 
V b (q ) =V b (q = 0) + q 



a + + a_ + 2c 

1 + a+ 
b + + &_ — a + — a_ 



(25) 
(26) 



Therefore, if we vary qo in the range —1/2 < go /eg < 1/2, 
we may tune it such that the equality in Eq. ( |2lj ) is satis- 
fied. In Fig. H we show the domain in the plane of Cq/C 
and qo/e where the SCB gap can be observed at zero 
temperature for the cases C x /C — and 1/4; the region 
enclosed by two thick solid lines represents the domain 
where the SCB gap is observed at zero temperature for 
C x /C — 0, and that enclosed by two thin solid lines for 
C x /C = 1/4. Fig. I indicates that for C /C > 0.1, the 
SCB region is readily accessible by adjusting the (gate) 
plate voltage V g . The region becomes wider as Co/C be- 
comes larger, but if Co/C is too big (> 10) the SCB gap 
shrinks, as does the usual Coulomb gap, and the current 
peak between them becomes too narrow to be observed. 
We thus find that the range of 0.1 < Co/C < 10 is prac- 
tically the range where the SCB gap can be observed, 
by adjusting the gate voltage V g . Fig. || also shows that 
the effect of the screening factor C x /C is to shrink the 
domain. That the domain shrinks with respect to the in- 
crease of C x /C may be understood by noting that, if the 
interaction between the top (2) and bottom (4) islands 
increases by the increase of C x /C, the stationary charge 
configurations are harder to be built up, so one needs 
less coupling between islands, by having bigger Co/C, to 
compensate the increase in the interaction between the 
islands. 

We can infer from Fig. || that, even when the back- 
ground charge qo/e of each island is induced, Co also 
increases as C x increases, as in the case of neutral back- 
ground charge (see Fig. 0). 
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FIG. 5. The domain in the plane of Co/C and qo/e where 
the SCB gap at zero temperature is observed, for C x /C = 
(the region enclosed by thick solid lines) and C x /C = 1/4 
(the region enclosed by thin solid lines). 
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FIG. 6. The temperature evolution of the SCB gap for 
Co/C = 1.0, C x /C = 0.25 and q /e = 0.3. The tempera- 
tures considered are T/T = (thick solid line), 0.0005 (thin 
solid line), 0.001 (dotted line), and 0.003 (dashed line). 



C. Effect of Fluctuations 

We have so far investigated the SCB phenomena when 
the temperature is zero and the quantum tunneling is 
suppressed. In the fluctuation-free case, the existence of 
the SCB gap results from the fact that the SCC with two 
trapped electrons is a stable configuration. The thermal 
and/or quantum fluctuations will, however, destabilize 
the SCC, so that the SCB gap will be transformed into 
a NDR region or entirely disappear, depending on the 
degree of the fluctuations. We will first discuss the effect 
of the temperature as follows. 

At finite temperature T, the trapped electrons in the 
SCC-2 have a probability T b to tunnel forward onto the 
island 3: 



1 



-AF h 



Re 2 1 - exp(AF 6 /fc B r) ' 



(27) 



The average time that it takes for the tunneling event to 
occur is (Atb) — l/2T b . (Tunneling out of both of the 
islands 2 and 4 contributes to the factor of two.) That 
is, after (Atb) on average, the SCC-2 is broken down by 
tunneling of one of trapped electrons onto island 3. The 
electron subsequently exits through the drain, provoking 
a series of tunneling events to take place before the SCC-2 
is restored, and the sequence of break-up and restoration 
of the SCC-2 is repeated. Let us denote (At a ) as the 
time that it takes for the SCC-2 to be restored after it is 
broken down and (AQ a ) as the amount of charges that 
are transfered through the array meanwhile. The current 
/ through the array is then given by 



I = (AQ a )/((At a ) + (At b )) 
= I a /(l + l/2T b (At a )), 



(28) 



where I a = (AQ a )/(At a ) is approximately the current 
that would flow through the array if the SCC-2 would 
not have been formed in the first place: i. e. the current 
which would fill the SCB gap smoothly (see Fig. ||). Since 
T b ~ AFb exp(— AF b /k B T) at low temperatures, 



I ~ I a exp(-AF b /k B T). 



(29) 



Therefore, for T < T c « AF b /k B , the SCB gap at zero 
temperature is transformed into a NDR region and for 
T > T c , the SCB gap entirely disappears and the current 
increases monotonically with respect to the bias voltage. 
See Fig. |. We estimate that T c /T = k B T c /(e 2 /C) w 
10 , which implies that for the junction capacitance C 
of the order of 1 aF, T c w 1 K. 

The second source of fluctuations is cotunneling. It 
is energetically favorable that the trapped electrons in 
the SCC-2 cotunnel onto the drain electrode. The co- 
tunneling rate j b is given by Eq. ( |l4| ) with AFi = AFb, 
AF 2 =a-V-(a+- 4c)/2 and AF = cV + (b+ + 36_)/2. 
Once the cotuneling event occurs, the same sequence of 
break-up and restoration of the SCC-2 takes place as in 
the case of the thermal fluctuation. The current / is 
then given by Eq. (28) with Tb replaced by j b . Since 
jb oc Rx/Rd, where Rd is the resistance of the junc- 
tion between island 3 and the drain, one roughly has 
I oc I a /Rd- Therefore, as in the case of the thermal fluc- 
tuation, there exist a certain Rd — R c above which the 
SCB gap is transformed into a NDR region and below 
which the I — V characteristics show monotonic increase. 
See Fig. 0. We estimate that R c w 10 5 fL 

It may be worthwhile to discuss on how cotunneling 
was incorporated in our calculation of the current in a lit- 
tle more detail. The current in the SCB gap was obtained 
by solving the master equation with the cotunneling pro- 
cess mentioned above included. As an alternative way, 
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FIG. 7. The evolution of the SCB gap with respect to the 
cotunneling strength at zero temperature. The array param- 
eters are the same as in Fig. ^| The resistances considered 
are Rd/R = 10 (thick solid line), 1 (thin solid line), 0.25 (dot- 
ted line), and 0.1 (dashed line), where Rd is the resistance 
of the junction between island 3 and the drain electrode and 
R = lMfi. 



we evaluated Eq. (|28|) with Tb replaced by jb and by in- 
dependently measuring I a through MC runs. Both of the 
results agreed very well. Using the cotunneling rate jb 
given by Eq. ( p^ ) in the two methods poses no difficulty 
in calculations in the SCB gap (V2 < V < Vb). Beyond 
the gap (V > Vb), however, 7^ diverges because of the oc- 
currence of sequential tunnelings in parallel. TheEefore 
we used the approximation by Jensen and Martinial3 to 
calculate the current beyond the gap and adjusted thus- 
obtained current curve so as for the current to be contin- 
uous at the threshold voltage Vb- This method of ours is 
far from being a rigorous one but may be justified because 
we are only interested in the effect of cotunneling in the 
SCB gap region for which we have treated the problem 
more carefully. 

Discussing both the thermal and the quantum fluctu- 
ations separately, we can now address the crossover tem- 
perature T x between the quantum-fluctuation dominant 
regime and the thermal-fluctuation dominant regime: for 
T < T x (the quantum regime), the currents are vir- 
tually unchanged with respect to the temperature but 
T > T x (the classical regime), the currents monotonically 
increase with increase of temperature. We may estimate 
T x by equating the normal transition rate Tb and the co- 
tunneling rate 7;, at a bias voltage V m for which we take 
the middle point of the SCB gap. Fig. || shows thus- 
estimated crossover temperatures versus the junction re- 
sistance Rd for the case of qo/e = 0.3, C x /C = 0.25, and 
Cq/C = 1.0. As is expected, as the junction resistance 
Rd increases, the crossover temperature T x decreases in 
the figure, which implies that the quantum fluctuations 
become less effective as Rd increases. 




RJR 



FIG. 8. The crossover temperature T x /T with respect to 
Rd/R- 



D. Effect of Disorders 

Let us now discuss how rigid the SCB gap is against 
possible imperfection of the array. Mainly there are three 
sources of disorder in real experiments: 1) the size of the 
islands may differ from one another, 2) the positions of 
the islands could be somewhat away from the desired 
ones, and 3) the stray charges may be induced on each 
island. Before we remark on the effect of the disorders 
in general, let us take a specific example to demonstrate 
the effect of the stray charges, among others. Note that 
stray charges on islands 2 and 4 may substantially affect 
the formation of the SCC-2 because the charges should 
be trapped there. When the background charge of island 
2 differs from that of other islands by Sq, the threshold 
voltage Vb changes by 



SV b = 



-Sq 
-Sq 



b--c 



a~ —< 
b+-c 



if Sq > 
if Sq < 0, 



and, regardless of the sign of Sq, V% changes by 
5V 2 = Sq 



I + 04 



(30) 



(31) 



The change in the threshold voltage Vt is immaterial in 
the discussion here. Eqs. (|3C|)-([3l|) show that the SCB 
gap shrinks linearly with respect to Sq. The SCB gap will 
be still seen if the change SVb — 8V2 < AVscb = Vb — V2, 
and we find that the range in Sq/e for which the SCB 
gap persists is typically around 0.1. But the SCB gap 
is much more vulnerable to negative stray charges than 
positive ones, which reflects the fact that a negatively 
charged electron is harder to get trapped if the island 
is more negatively charged. If we can adjust the overall 
background charges on each island, in the direction to 
compensate the stray charge, we may be able to re-enter 
the region where the SCB gap is seen. 
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We found it a challenging task to treat the problem 
analytically when the disorders are present in general, so 
we resorted to numerical simulations with the following 
simple model. The first two sources of imperfection men- 
tioned above may be treated by relaxing the condition of 
having identical islands and allowing the mutual capaci- 
tances to have 'random' contributions to some degree as 
follows: 

Qj^CijO-aCu) for i^j, (32) 

where CV,- refers to the capacitance between island i and 
j (including the source, the ground, and the drain elec- 
trodes), a is a constant adjusting the magnitude of the 
randomness and £'s are random numbers between -1/2 
and 1/2. Note that in the above equations, we have 
specifically put the subscripts for £ to note that the ran- 
dom numbers are differently assigned for different islands 
and different pairs of islands. Similarly, for the third 
source of imperfection, we may introduce random contri- 
butions to q of each island. Our numerical simulations 
show that for up to 10 % of random contributions the 
SCB phenomena is still observed, although the SCB gap 
position or the NDR region is shifted and/or its width is 
changed. 

IV. CONCLUSION 

In conclusion, we have shown that the secondary 
Coulomb blockade gap of the four-island array at zero 
temperature is originated from the property of the array 
that it traps up to two electrons inside, due to its unique 
topology. The detailed tunneling mechanism which leads 
to the formation and destruction of the stationary charge 
configuration with two trapped electrons inside, which 
corresponds to the additional gap, has been analyzed. 
We have also considered the effect of the thermal and 
quantum fluctuations on the gap, and found that the 
gap transforms into a NDR region upon introduction of 
the fluctuations. The NDR feature suggests a potential 
application of the ring-shaped array as a diode showing 
the NDR behavior. We have also studied the effect of the 
uniform background charges on each island, which can be 
realized in real experiments by attaching a ground plate 
near the array and applying a voltage, and found that 
the secondary Coulomb blockade gap is more readily ac- 
cessible by adjusting the voltage. Based on the detailed 
analysis that we have carried out in this study, the physi- 
cal properties of the multiple Coulomb gaps of larger size 
array may be similarly understood. 

The authors thank Gwang-Hee Kim for fruitful discus- 
sions on the cotunneling effect and the crossover temper- 
ature. This work has been supported by the Ministry of 
Information and Communications of Korea. 
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